Single Spin Dynamics and Decoherence in a Quantum Dot via Charge Transport 
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We investigate the spin dynamics of a quantum dot with a spin-i ground state in the Coulomb 
blockade regime and in the presence of a magnetic rf field leading to electron spin resonances 
(ESR). We show that by coupling the dot to leads, spin properties on the dot can be accessed via 
the charge current in the stationary and non- stationary limit. We present a microscopic derivation 
of the current and the master equation of the dot using superoperators, including contributions to 
decoherence and energy shifts due to the tunnel coupling. We give a detailed analysis of sequential 
and co-tunneling currents, for linearly and circularly oscillating ESR fields, applied in cw and pulsed 
mode. We show that the sequential tunneling current exhibits a spin satellite peak whose linewidth 
gives a lower bound on the decoherence time T2 of the spin-| state on the dot. Similarly, the spin 
decoherence can be accessed also in the cotunneling regime via ESR induced spin flips. We show 
that the conductance ratio of the spin satellite peak and the conventional peak due to sequential 
tunneling saturates at the universal conductance ratio of 0.71 for strong ESR fields. We describe 
a double-dot setup which generates spin dependent tunneling and acts as a current pump (at zero 
bias), and as a spin inverter which inverts the spin-polarization of the current, even in a homogeneous 
magnetic field. We show that Rabi oscillations of the dot-spin induce coherent oscillations in the 
time-dependent current. These oscillations are observable in the time-averaged current as function 
of ESR pulse-duration, and they allow one to access the spin coherence directly in the time domain. 
We analyze the measurement and read-out process of the dot-spin via currents in spin-polarized 
leads and identify measurement time and efficiency by calculating the counting statistics, noise, and 
the Fano factor. 
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I. INTRODUCTION 



The coherent control and manipulation of the elec- 
tron spin has become the focus of an increasing num- 
ber of experiments. 1-8 From measurements it has be- 
come evident that the phase coherence of electron spins 
in semiconductors can be robust over-, unusually long 
times, exceeding 100's of nanoseconds. a Thus, spins of 
electrons are suitable candidates for applications in the 
field of spintronics, in particular for quantum information 
processing. 9-18 This has made it desirable to understand 
in more detail the coherent behavior of single electron 
spins which are confined to nanostructures such as quan- 
tum dots, molecules, or atoms, and to point to ways of 
how to access the coherence time of a single spin experi- 
mentally. It is the goal of this work to address this issue 
and to propose and analyze transport scenarios involv- 
ing a quantum dot attached to leads and with a spin-1/2 
ground state. 

Wc first remind ourselves of some basic notions in spin 
dynamics. When the electron spin is exposed to an exter- 
nal magnetic field, this leads to a Zeeman splitting, and 
the spin dynamics is described by the standard Bloch 
equations.E3 These are characterized by two time scales: 
the (longitudinal) relaxation time T\ and the decoherence 
time T2 (transverse relaxation) . The spin relaxation time 
T\ describes the lifetime of an excited spin state, aligned 
along the external field, and is classical in the sense that 
its definition does not involve the concept of quantum su- 
perpositions. Such a T\ time of a spin in a single quantum 
dot was measured recently via transport and was shown 
to be longer than a few microseconds ,0 in agreement with 



calculations .0 On the other hand, the spin decoherence 
time T2 gives the time over which a superposition of op- 
posite spin states of a single electron remains coherent. 
Thus, coherent manipulations of electron spins, e.g., gate 
operations for quantum computation, must be performed, 
faster than T%. We note that quite generally T2 < T\x3 
Thus, from the sole knowledge of T\ no lower bound for 
T2 follows. It is thus of fundamental interest to investi- 
gate possibilities of how to gain access to the decoherence 
time T2 for a single spin confined to a quantum dot. 

The loss of phase coherence of many but independent 
spins is described by the dephasing timed T 2 *, where in- 
homogeneities in the Zeeman terms lead to a further sup- 
pression of phase coherence for the ensemble but not nec- 
essarily for an individual spin, thus T 2 * < T2. In recent 
experiments, T 2 * was measured in bulk GaAs by using 
ultrafast time-resolved optical methods, yielding values 
for T| exceeding 100 ns.d 

However, the measurement of the decoherence time T2 
for a single spin has — to our knowledge — not been re- 
ported yet (although it is expected to be within experi- 
mental reach given the known single-photon sensitivity). 
A first step into this direction are spin echo measure- 
ments on an ensemble of spins, where dephasing due 
to inhomogeneities of the magnetic field is eliminated. 
Indeed, such measurements being performed more than 
thirty yeacs ago on P donors in Si, reported T2 times up 
to 500 fisE3 However, it appears desirable to have a more 
direct method for single spin measurements. To achieve 
this via direct coupling to the magnetic moment of the 
spin is rather challenging due to the extremely small mag- 
netic moment, although it is-believed to be within reach 
using cantilever techniques E3 Here we concentrate on a 
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further approach based on transport measurements. The 
key idea is to exploit the Pauli principle which connects 
spin and charge of the electron so intimately that all spin 
properties can be accessed via charge and charge, cur- 
rents, especially in the Coulomb blockade regimecj of a 
quantum dot attached to leads. Indeed, concrete scenar- 
ios based on such a sjmjJ^tCjharge conversion have been 
proposed in the past,LflljO'E£tE£l and it is our goal here to 
further elaborate on these concepts, and to report on a 
variety of new results we have obtained. 

There are two classes of spin decoherence contributions 
we have to distinguish in the following. First, rare tun- 
neling events of electrons onto and off the dot change the 
spin state of the dot and in this way contribute to the 
decoherence of the dot spin. We account for this deco- 
herence microscopically in terms of a tunneling Hamilto- 
nian. Second, there are intrinsic decoherence contribu- 
tions from processes which persist even if the dot is com- 
pletely isolated from the leads. This decoherence is taken 
into account phenomenologically in the master equation 
developed in this work, with an intrinsic decoherence rate 
T^ 1 . The goal then is to show that this T2 time can be 
extracted via current measurements, regardless of the mi- 
croscopic processes leading to T2. Such a phenomenolog- 
ical approach to intrinsic decoherence makes the purpose 
of our considerations clearer and is applicable to different 
types of decoherence mechanisms e.g. based on hyperfine 
and spin-orbit couplings. The microscopic study of such 
intrinsic decoherence, being an important subject in its 
own right, is not addressed in the present work. 

The outline of this paper is as follows. In Sec. II, 
we define the system of interest, a quantum dot with 
spin- 1/2 ground state in the Coulomb blockade regime 
tunnel-coupled to leads, and in the presence of an elec- 
tron spin resonance (ESR) field. We derive the (gen- 
eralized) master equation for the low energy dot states 
in the sequential and co-tunneling regime by evaluating 
the tunnel coupling to the leads microscopically in order 
to obtain tunneling rates, decoherence rates, and energy 
(Stark) shifts. For this we need to include diagonal and 
off-diagonal matrix elements of the reduced density op- 
erator. The stationary current through the dot and its 
dependence on the ESR field is discussed in Sec. [II. Wc 



find a spin satellite peak in the sequential tunneling cur- 
rent, whose linewidth as function of the ESR frequency 
gives a lower bound for the T2 time. Thus, via the sta- 
tionary current, the T2 time can be accessed in a regime 
that is experimentally accessible, as will be demonstrated 
by concrete numerical examples. We show that the ra- 
tio of this satellite peak and the main peak saturates at 
a universal conductance ratio for strong ESR fields. In 
Sec. IV, we extend our results to the even-to-odd tran- 
sition, i.e., for the case where there is (on average) one 
electron less on the dot. In Sec. we explain a mech- 
anism for a spin inverter device which inverts the spin- 
polarization of the current passing through two dots cou- 
pled in series in the presence of a homogeneous magnetic 
field. In Sec. VI, we discuss how spin-dependent tunnel- 



ing can be used to pump a current through a system in 
the absence of a bias, w here the ESR field provides the 
required energy. In Sec. VII, we consider rotating ESR 



fields which allows us to obtain the exact time evolution 
of the dot-states and their decay rates. In Sec. 



VIII, the 



cotunneling current through the quantum dot away from 
the sequential tunneling peak is discussed. We show that 
the T2 time can also be accessed in this regime. Invoking 
spin-polarized leads, a read-out procedure for the dot- 
spin is proposed and analyzed in Sec. IX, where count- 
ing statistics, noise, and the Fano factor are calculated, 
which allow us then to estimate the measurement time. 
In Sec. |x|, we discuss coherent Rabi oscillations of the dot 
spin and their occurrence in the time-dependent current. 
In Sec. XI , we show that Rabi oscillations can also be ob- 
served in the time-av erag ed current if pulsed ESR fields 
are applied. In Sec. XII , we point out that our results 
also appl y to STM devices, and we finally conclude in 



Sec. XIII 



II. QUANTUM DOT IN ESR FIELD 

A. Model Hamiltonian 

We consider a quantum dot in the Coulomb blockade 
rcgime,Ea which has a spin-i ground state. The dot is 
assumed to be tunnel-coupled to two Fermi-liquid leads 
Z = 1, 2, at chemical potentials fii. We start from the full 
Hamiltonian 



H = H\ 



lead 



#dot + #ESR.(*) + #T> 



(i) 



which describes leads, dot, ESR field, and the tunnel cou- 
pling between leads and dot, respectively. For the leads 
we take #i cad = £ /fc(T e lk cj k(7 c lkl7 , where c\ k<7 creates an 
electron in lead I with orbital state k, spin a, and energy 
eik- We describe the coupling with the standard tunnel 
Hamiltonian 



H T = 



Ipka 



h.C. 



(2) 



with tunneling amplitude tf p , and where d P<7 creates an 
electron on the dot in orbital state p. In Eq. ([!]), H& ot 
is time-independent and includes charging and interac- 
tion energies of the electrons on the dot and coupling to 
a static magnetic field B z in z direction. The dot-spin 
is coupled to a magnetic ESR field, B x (t) — £?"cos(wi), 
linearly oscillating in x direction with frequency w, thus 
Hesr = —^g^BB x (t) a x . Such an oscillating field pro- 
duces Rabi spin-flips when its frequency is tuned to reso- 
nance, uj = A z , as shown below. Then, the total Zecman 
coupling of the dot-spin is 



-^ B B(i) • er 



1 



1 



A^cos (0^)0-3;, (3) 



with electron g factor g, Bohr magneton //b, and Pauli 
matrices <r,-. We have defined = 5/13-6°, and the 
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FIG. 1: Quantum dot coupled to (unpolarized) leads I = 
1, 2 with chemical potentials /i;. The sequential tunneling 
regime Es > fii > Es — A z > fi2 (for i? T = 0) shown here 
corresponds to t he sat ellit e pea k in the sequential tunneling 

and Figs, u and H. Here, Es 



III A and III B 



current, cf. Sec. 
(Et + ) are the singlet (triplet) levels and the Zeeman splitting 
is A z = gr/iB-Bz > kT. (a) If the dot is initially in the spin 
ground state | |) , sequential tunneling is blocked by energy 
conservation, (b) If the dot-spin is excited by an ESR field 
(Rabi flip), spin up electrons can tunnel from lead 1 onto the 
dot, forming a singlet. Then, spin up or down electrons can 
tunnel into lead 2. 



Zeeman splitting A z = <?/iB-3 2 . Ideally, we assume that 
the Zeeman splitting of the leads A' 2 cads is different from 
A z , and A' 2 cads -C £f, where £p is the Fermi energy, such 
that the effects of the fields B z and B x (t) on the leads are 
negligible (see below). Such a situation cap be achieved 
by using materials of different g factorsEl and/or with 
local magnetic fields (B x or B z ). 

We arepiieglecting photon assisted tunneling (PAT) 



processeSjCilcZI in which oscillating electric potentials of 
the leads provide additional energy to electrons tunnel- 
ing onto the dot. We note that PAT contributions to the 
current can be distinguished from ESR effects since the 
former contributions do not show resonant behavior as a 
function of B z and/or u>, and they lead to several satel- 
lite peaks instead of one as for ESR effects (see below). 
Further, if one avoids electrical rf components parallel 
to the current, i.e., along the axis lead-dot-lead, no po- 
tential oscillations are produced, and thus PAT effects 
are excluded. Finally, electric rf fields can be avoided 
altogether, using a setup as in Ref. [2^. There, the oscil- 
lating current induced in a superconducting wire (via an 
rf source) generates only a magnetic rf component in the 
near-field region,E3 with an the electric component that 
is negligibly small for u> <C u> p , where w p is the plasma 
frequency. 



B. Dot Spectrum and Energetics 

The electronic states of the quantum dot can be as- 
sumed as follows. For an odd number N of electrons on 
a dot with antiferromagnetic filling, the dot has a spin- 
\ ground state. The topmost (excess) electron can be 



either in the spin ground state, IT), {p z eigenstate) or 
in the excited state, | J.) (see Fig. |l|). This assumption is 
automatically satisfied if i\T = 1. Otherwise, to obtain an- 



tiferromagnetic filling, Hund's rule must not apply. This 
can be achieved by breaking the orbital degeneracy on 
the dot, e.g., by using asymmetrically shaped dots or an 
appropriate magnetic field -B 2 .Ej For an additional elec- 
tron on the dot, we assume for N + 1 the ground state 
to be the singlet \S) = (|||) - \IW/V%, i-e., the triplet 
state \T + ) = 1 1|) has higher energy, which again can be 
achieved by tuning i? z .Ej The energy E m of the dot, in- 
cluding charging energy, is defined by Hd ot \m) = E m \m). 

We shall give a brief overview of the energetics involved 
in tunneling through quantum dots in the Coulomb 
blockade regimes and in the presence of the Zeeman 
splitting and an ESR field. For simplicity, we assume 
that there is no electron-electron interaction on the dot 
apart from the classical charging effect. (Our work is not 
restricted to such an assumption, since we only require 
a spin-i ground state and a large enough singlet-triplet 
spacing on the dot.) The total ground state energy of a 
dot with antiferromagnetic filling is 



U(N) 



N 

fc=l 



E, 



N 



(4) 



for N electrons on the dot. Here, the single particle en- 
ergy of the fcth electron, e a k = eu + (— l) fc A 2 /2, contains 
orbital and Zeeman energy contributions. The charging 
energy is E^ = (Ne — Qg) 2 /2C, with gate charge Qc, 
and dot capacitance C. It is convenient to define the 
chemical potential of the dot, [idot(N + 1) = U (N + 1) — 
U(N), which is the energy required for an electron of lead 
I to tunnel onto the dot, which contains N electrons ini-, 
tially, i.e., tunneling onto the dot occurs for /x; > /idot-" 
In the Coulomb blockade regime, kT <C e 2 /C (k: Boltz- 
mann constant), no sequential tunneling current flows 
through the dot if the chemical potentials of dot and leads 
are such that Hdot(N) < [x%, [12 < ^dot{N + 1). However, 
in the sequential tunneling regime [i\ > /J-dot{N+l) > (12, 
single electrons tunnel from lead 1 onto the dot and then 
on into lead 2, producing a sequential tunneling current. 

In the presence of an ESR field, these concepts must 
be extended. Excitations of the dot states must be taken 
into account, since now the energy of the dot changes in 
time due to B x (t). A full analytical description of the 
current flow is derived in the following sections based on 
a time dependent master equation. Here, we just intend 
to give a qualitative picture to provide some intuition for 
the underlying physical mechanism (it will not be needed 
later on) . We define a time-dependent chemical potential 
of the dot, given as the energy required to add an electron 
at time t. We consider the two chemical potentials [i^ot 
for initial spin-i dot-state \a), i.e., Agf = n\ot(N + 1) = 

Es —E\, which 
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Es — A z , respectively, 



Es -E h and A sl = fi dot {N - 
simplify to Asi — Es, and Asi 
for Ei = 0. Note that the ^ ot is lowered if the dot is 
excited into state ||), since the Zeeman energy A 2 has 
already been provided by a Rabi spin flip due to the 
ESR field. Therefore, we can identify the regime Agj > 
fi± > Asi > A^j where a sequential tunneling current 
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will flow through the dot only after exciting the dot-spin 
by a spin flip (see Fig. [I]). In other words, the dot can be 
opened and closed via the ESR field, which thus allows to 
modulate the current. This (dynamical) dependence of 
the current on the dot-spin can be exploited to measure 
the Ti time and the Rabi oscillations of the dot-spin,E£l 
as we will explain in detail in the following. 



C. Systematic Treatment of Sequential Tunneling 

The electronic states on a quantum dot interact with 
their environment (heat bath), in particular with the 
Fermi leads, which provide and take up electrons. The 
state of the combined system, dot and environment, is 
given by the full density matrix pit). The states of in- 
terest are the electronic states on the dot, described by 
the reduced density matrix of the dot, pd = Ttb p. Here, 
Trs is the trace taken over the leads (environment), av- 
eraging over the (unobserved) degrees of freedom of the 
environment. The diagonal elements p n = (n\pu\n) of 
the density matrix of the dot describe the occupation 
probabilities of the dot levels, with H^otW) = E n \n). 
The off-diagonal elements p nm = (n\pY>\m) — p* nn de- 
scribe the coherence and the phase of superpositions of 
dot-states. 

The tunnel coupling Ht between leads and dot is 
switched on at t = 0. Prior to this, the dot and leads 
are assumed to be uncorrelated such that the full initial 
density matrix factorizes as p(0) = pd(0)p^, where p^ is 
the density matrix of the leads in thermal equilibrium at 
/ii.2, and at temperature T. Next we derive the master 
equation for the reduced density matrix pd by making 
use of the superoperator formalismEj In the following, 
we set h = 1. Starting from the von Neumann equation 
p = —i [H, p] for thc-jjull density matrix, and using stan- 
dard manipulations ,E3 one finds the time evolution of the 
reduced density matrix 



Po{t) = -i [H dot + H ESR (t), pn(t)] 



dt'M(t, t')pn(t') 



(5) 



M (t, t') 



Tr B L T Te 



LtpI, (6) 



with time-ordering T and the Liouville operators defined 
by L(t)X = [H(t), X], L T X = [H T , X], and equiva- 
lent^ for L dot , L lcad , and L ESR (i). The projectors are 
defined as Q = 1-P, and PX = p^Tr B X. The kernel M 
[Eq. (^j|)] is a superoperator describing processes involving 
tunneling of electrons to and from the leads. We consider 
here only sequential tunneling processes, and re fer f or a 
discussion of cotunncling contributions to Sees. [IF and 
VIII. Thus, we work in Born approximation by retaining 
only the terms in lowest order of Lt, i.e., we replace L by 
Lq = L — Lt in Eq. (^|). For further evaluation of M, it is 
self-consistent (see below) to neglect the effect of the ESR 
field, L ESK (t), i.e., we replace L by L dot + Li cad in M. 



This removes explicitly the time dependence of M, mak- 
ing it time translation invariant, M(t, t') — M[t — t'). 
We find that M(r) decays on a time scale r c ~ 
i.e., the correlations induced in the leads by Ht decay 
rapidly. Since this decay is typically much faster than 
the Rabi flips produced by the ESR field, r c <§; 1/A X , 
we may indeed neglect the contribution of Lesr(£) to 
M . With these approximations, Eq. (|J) becomes in the 
interaction picture 

p\{t) = -iI% SSR (t)fk(t) - fdTM\T)p\,{t-T).{7) 

Jo 

The rapid decay of M(r) also justifies the Markovian 
assumption that the system has no memory about its 
past, i.e., that f>T>{t) depends only on pj)(t) and not on 
Pn{t — t). This approximation is performed in the in- 
teraction picture, to keep track of the dynamical phase 
of the off-diagonal elements of pn- Systematically we 
proceed as follows. Since the integrand in Eq. (0) only 
contributes for small t, we may expand the integrand 
in r, M(ry D (i - r) = M{r)[p\{t) - rp\{t) + 0(r 2 ]J. 
We then replace p^it) in the integrand by using Eq. (|7) 
iteratively. However, since M(t) - 0(L T 2 ), we can ne- 
glect the part of Pq(<) which is 0(Lt 2 ), since it corre- 
sponds to a higher order term in our Born approxima- 
tion. The remaining part of p\y{t) results from LesRj 
which can also be disregarded since, in the integrand, 
the ESR field only acts on the time scale r c ^ 1/A K . 
We then extend the upper integration limit in Eq. (j?]) 
to oo, with negligible contributions due to the decay of 
M(t). Therefore, the second term in Eq. (^|) becomes 

— { f^dr M 1 (r) } pj-, (t) . Next, we evaluate the matrix el- 
ements M bc \ nm = {b\ f M \n){m\) \c) explicitly in the in- 
teraction picture, which yieldso 

- dr mI~(t) = S bc 6 mn ( W cn ~ S bn V W kn 

-(1 - 5nm) $ 

k 

with the rates W (see below) and energy shifts Se nm 
(Stark shifts). These shifts are small; e.g., the one be- 
tween 1 1) and 1 1) is given by 



1 f°° 



7, 



7; 



e— A 



5'1 



£ - A 



s | 



, (9) 



and similarly for Sesi and Ses^- For |p; — A^o-I > kT, 
the energy shift becomes 



Se 



LT 




(1, 



A 



2l 

2k 



log 



(1, 



A 



si 



T 

which, for 7^ = 



(10) 



7; 



reduces 



to 6e iT « 

E/(7;/27r)log [J/xj - Asxl/|p; - A ST |], and thus to 
a small correction |(5ejj | < 7log(A^/fcT), for Ap < A z . 
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The sequential tunneling rates in Eq. (||) are 

i 

W lS = J2 W ls> W ls = 7/ [1 - M&si)] , (12) 



with the Fermi function fi(A Si ) = [l + e {A ^-^/ kT ] 1 
of lead I. The rates Wgf, Wfs, W S p and are de- 
fined analogously as functions of 77 and fi(As^). The 
transition rates, 



7 ; = 2^*; 



n 2 



7; 



27T^| f 



11^ 



! I 2 



(13) 



consist of (possibly) spin-dependent density of states 1^ 
at the Fermi energy and tunneling amplitude t, . (Spin- 



dependent density of states are considered in Sec. IX for 
spin read-out.) For later convenience, we define for a =| 

, i 



7 CT = (7i CT +72)/2, 7 =(7 T +7 X )/2 



D. Master Equation 



(14) 



So far we have considered only coupling to an envi- 
ronment consisting of Fermi leads. However, the elec- 
tronic dot states are affected also by intrinsic degrees 
of freedom such as hyperfine coupling, spin-orbit inter- 
action, or spin-phonon coupling, which lead to intrin- 
sic spin relaxation and decoherence. Treating such cou- 
plings microscopically is beyond the present scope (see 
e.g., Ref. |2(]). Thus, we treat these couplings phenomeno- 
logically by introducing corresponding rates in the mas- 
ter equation. First, the spin relaxation rates Wfi and 
W11 describe processes in which the dot-spin is flipped. 
We can assume Wfi >• W^, for A z > kT (consistent 
with detailed balance, W^/W^ = e A ^ feT ). These re- 
laxation processes correspond to the phcn omcnological 
rate 1/Ti = W n + W^, see also Sec. |ll| Second, the 
rate 1 /T2 describes the intrinsic decoherence of the spin 
on the dot, which is present even in the absence of cou- 
pling to the leads. This type of decoherence destroys the 
information about the relative phase in a superposition 
of |t) and I J,), without changing the populations of the 
opposite spin states. Formally, this leads to a decay of 
the off-diagonal matrix element p it . Including the deco- 
herence contribution of Ht [Eqs. (||), the total spin 
decoherence rate is 



W S ] + W sl 



1 



(15) 



i.e., electrons tunneling onto the dot further destroy spin 
coherence on the dot (see Sec. [IE for an interpretation). 

With the above results, we obtain from Eq. (||) the 
master equation of the dot, 



Pi 



Pi 



-(Wit + Ws]) p T + W n pi + W 1S Ps 
-A x cos (cjt) Im[p iT ] , (16) 
VF iT p T - (W n + W Si ) Pi + W IS p s 
+A a cos(a;i)Ini[p| T ], (17) 
IFstPt + w siPi ~ {Wis + W ls )p s , (18) 



Ps 

PIT = - iA zPl1 + i^-cos(uit)(p 1 - Pi ) - %p iT ,(19) 



Ps] — —iA. S ]Ps\ — Vs]Psii 
Psi = -iAsiPsi - Vsi psi- 



(20) 
(21) 



Here, the time evolution of the matrix elements p nm = 
(n\pn\m) of the density matrix of the dot is described 
for the states \n) = II), \S), e.g., for the diagonal 
element we write pj = (f |pd| T), and for the off-diagonal 
element, psi — (S\pu\1), etc. The rate W mn describes 
transitions from state \n) to |m). Equations (16)-jT^) 



are rate equations with gain and loss terms, up to the 
contributions from the ESR field. Then, the population 
of, say, state 1 1), is changed by dp^ after time dt by the 
following contributions [Eq. (jL6|)]. The population pf is 
increased when the dot is previously in state IS") (with 
probability ps), and a spin J, electron tunnels out of the 
dot with probability W^s dt. However, the population pj 
is decreased when the system was already in state |t), 
and a spin J, electron tunnels onto the dot with proba- 
bility Ws\ dt. The spin flip rates, Wfi and Wj,j, enter 
Eq. (|l^) analogously. In the absence of an ESR field, 
the off-diagonal elements [Eqs. (|l9|)-(pT|)] of the density 
matrix decouple from the diagonal ones and decay with 
the decoherence rates V nm — V mn . 

In the presence of an ESR field, the diagonal [Eqs. ( fl6| ) 
and ([IT])] and the off-diagonal [Eq. (|l9|)] matrix elements 
become coupled by the term proportional to A^. This 
coupling of populations (pf and p[) and coherence (pxt) 
shows the coherent nature of Rabi spin-flips and makes it 
apparent that we are studying a resonant process, which 
requires that we take H^sr fully into account. 

The current I2 = (dq/dt) from the dot into lead 2 is 
defined by the number of charges dq that accumulate in 
lead 2 after time dt. With probability ps, the dot is 
in state \S) and a charge will tunnel into lead 2 with 
probability (W? s + Wf s ) dt. However, if the dot is in 
state |t) or | j), a charge may tunnel from lead 2 onto 
the dot, reducing the number of charges in lead 2. Thus, 
in total we obtain for the current in lead 2 

h = e (Wf s + Wf s ) PS -eW 2 s1 pi-e W§ lPl . (22) 

The current in lead 1 , 1\ , is obtained analogously and is 
given by Eq. (|2^) after changin g sig n and replacing the 
index 2 by 1. We show in Sec. HI that l\ = I2 in the 



stationary limit, due to charge conservation. 

Finally we note that Eqs. (^) and (|l|), which de- 
scribe a superposition of an odd and an even number of 
electrons on the dot, decouple from Eqs. (16)-(|l9|) and 
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are thus not of relevance for our considerations. Fur- 
ther, since the coupling to the leads is switched on only 
at t = 0, initially the number of particles on the dot is 
well defined. Therefore ps] and psi vanish at t = and 
at all later times, as seen from Eqs. ( p0| ) and (|2l|). In 
particular, no superposition of a state with an even and 
a state with an odd number of electrons on the dot is 
produced by the coupling to the leads, since this would 
require a coherent superposition of corresponding states 
in the leads, however, for times larger than t c (which is 
typically the case), we can safely neglect any coherence 
in the Fermi liquid leads. 



E. Decoherence and Measurement Process 



F. Cotunneling Contribution to the Sequential 
Tunneling Regime 

We work in the sequential tunneling regime, defined 
by pi > Ag± > p2- One can see that hisbp, order — 
cotunneling — contributions can be neglect edE3EJ for 7; < 
A z , kT, the regime of interest here. Most importantly, 
the cotunneling con tributions to Vjj are of the order 
J?/A z (see Sec. VTl|), i-e., they are suppressed compared 
to the sequential tunneling contributions by a factor of 
ji/A z (« 5 x 10~ 5 for the parameters of Fig. [|). For- 
mally, the cotunneling contributions to the master equa- 
tion can be absorbed into T\ and T^. For a discussion of 
cotunneling curren ts aw ay from the sequential tunneling 
resonance see Sec. VIII. 



We elucidate the connection between spin decoherence 
and measurement, first in the absence of leads and ESR 
field. We consider a coherent superposition a\ f) + fi\ [) 
as the initial state of the dot. This pure state cor- 
responds to the reduced density matrix p|(0) = \a\ 2 , 
= I/3| 2 j and pj,|(0) = af3* , and the master equa- 
tion contains only the rates W^i, W|j,, and Vjj = I/T2. 
The off-diagonal terms pjj = pt,, decay with the de- 



coherence time T2, Pi^{t) 



- „-t/T 2 -itA z 



P1t(°)> wm l c 



the diagonal terms (occupation probabilities) decay with 
the spin relaxation time Ti = (Wfj + VF^) -1 , Pi(t) — 
<°i q + e~*/ Tl [/?x(0) — P? q ]j toward their stationary value 
PT = W u /(W n + W lr ), and p T = 1 - Pl . In total, for 
T2 < T\ , we can picture the decay of pr> as 



\a\ 2 a/3* 
a*(3 \(i\ 2 







CO 

p* 





eq 

Pi 



(23) 

i.e., the off-diagonal terms vanish first on the timescale 
T2, and then the diagonal ones equilibrate on the 
timescale T\. 

when electrons tunnel onto the 



As shown in Sec. II C 



dot, the decoherence rate Vj.| [Eq. (15)] and thus the de- 
cay of the off-diagonal elements is increased further. We 
note now the formal equivalence to the quantum mea- 
surement process (in the u z basis), where the dot-spin is 
projected onto | \) or | j), and thus the off-diagonal ma- 
trix elements vanish. This projection can be understood 
as a decoherence process. Conversely, we can consider 
the decoherence due to tunneling as a measurement per- 
formed by the tunneling electrons. We note that this pro- 
cess is a weak measurement in the following sense. The 
electrons in the leads attempt to tunnel on the dot, but 
only with small probability cx Wsa are these attempts 
successful. Thus, the current /, which carries away the 
information of the dot state to the observer, is formed 
by these successful electrons, while the unsuccessful elec- 
trons are not detected. Another way to say this is that a 
given electron from the lead has only a small probability 
oc Wsa to "measure" (i.e., decohere) the dot state. 



III. STATIONARY CURRENT 

We now consider the stationary current / in the pres- 
ence of a continuous wave (cw) ESR field. Therefore we 
calculate the stationary solution p(t — > 00) of the mas- 
ter equation [Eqs. (|l6|)-(pil|)l- We will apply the rotating 
wave approximation (RWA) ,Lj where only the leading fre- 
quency contributions of H^sr are retained. Higher order 
contributions would include simultaneous absorption of 
two photons and emission of another photon. In lowest 
order, only single photons can be absorbed or emitted, 
producing a spin flip on the dot. To perform this ap- 



proximation, we write A^cos (uji) = hA x (e 1 ' 



i.e., we decompose the linearly oscillating magnetic field 
into a superposition of a clockwise and an anti-clockwise 
rotating field. Integrating Eqs. (|l^), (|l^), and (|l^), one 
finds that for oj « A 2 , the anti-clockwise rotating field 
leads to rapidly oscillating terms in the integrands, which 
nearly average to zero. Therefore, we retain only the 
clockwise rotating field, which is given by the term pro- 
portional to e wt (see also Sec. VII). Note that since only 
one field component contributes, the field amplitude is 
halved. This leads to the period Tsi of one Rabi oscilla- 
tion, 



To 



47T 

a!" 



(24) 



The RWA is valid for A x , Vj.t , \ A z -u\ < u, see e.g., Ref. 
|35| , and is well justified for the parameters considered 
here. In the stationary case and using the RWA, the 
dependence of p-\ and p^ [Eqs. (|ll]) and (pT|)] on p^ is 
eliminated, leading to the effective spin-flip rate 



A., 
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A z y 



nv 



(25) 



which is a Lorentzian as function of lu with maximum 
W™ ax = A x 2 /8Vtf at resonance uj = A z . 

Now it is straightforward to find the stationary solution 
of the effective rate equations for /5f , p^ and ps, 

PT = ri[WisWsi + (W n +W„)(Wis + W ls j\, (26) 
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Pi = V [W iS W S] + (W u + W u ) (W 1S + W ls )] , (27) 

Ps = v [WsiWsi + Wsi(W n 

+W SI (W U +W U )], (28) 

where the normalization factor r\ is such that ^ p n — 1. 
We see from Eqs. (p6|)-(|28|) that the effective spin flip 
rates are W n + W^, and VK iT + W u , i.e., the ESR field 
flips up and down spin with equal rate W u . 

We can now calculate the spin-"f polarized current in 
lead 2, 4 = eWf s ps -eW^pi [cf. Eq. ©]. The re- 
sult is displayed in Eq. (Al) in the Appendix. The spin-j 



polarized current, Jj, is obtained from Eq. (Al) by in 



terchanging | with j in the numerator (the denominator 
remains unaffected by such an interchange). The currents 

t I T I 

in lead 1, I{ are obtained from the formulas for Jg by 
changing sign and interchanging indices 1 with 2. Note 
that generally l\ ^ l\ , since the ESR field generates spin 
flips on the dot, and thus the spin on the dot is not a con- 
served quantity. However, the stationary charge current 
II = If is the same in both leads, I = I\ = I2, due 
to charge conservation. 



A. Spin Satellite Peak 

In this subsection we discuss the stationary current / 
through the dot, in particular, its behavior as function of 
p = (/ii + /i 2 )/2, or, equivalently, as function of the gate 
voltage V g . We will see that an additional sequential 
tunneling peak (satellite peak) will appear due to the 
ESR field. Before explicit evaluation of the current, we 
briefly describe this situation in qualitative terms. We 
assume a large Zeeman splitting, A z > Ap, kT, with 
applied bias Ap = p\ — pi > 0. If the potentials are 
such that p\ > Asf > P2, i.e., the chemical potential 
of the dot (relative to the ground state ||)) is between 
the chemical potentials of the leads, the state on the dot 
changes between | f ) and \S) due to sequential tunneling 
events, leading to the standard sequential tunneling peak 
in I(p) at p w Asf . 

However, we also have to consider the regime Agj > 
pi > A s 1 > P2, as shown in Fig. Q. Without ESR field, 
the dot relaxes into its ground state 1 1) (since <C 
Wfj.), and the sequential tunneling current through the 
dot is blocked since the chemical potential Agf of the 
dot is higher than those of the leads. However, if an ESR 
field generates Rabi spin- flips (on the dot only), the cur- 
rent flows through the dot involving the state | J,), since 
Asi is lower than p±. Therefore, a sequential tunneling 
current appears also for gate voltages V g corresponding 
to As i, i.e., I(p) exhibits a spin satellite peak due to the 
ESR field at p w Ag^. This new peak is shifted away 
from the main peak by A 2 (Fig. ||). The presence of 
such a satellite peak and its sensitivity to changes in B z 
allows identification of spin effects.cJ Further, we note 
that via the position of the peak in I(u>), I(B Z ), or I(p), 
the Zeeman splitting and also the g factor of a single dot 




FIG. 2: The stationary current / [Eq. jAl[ )] vs. p = (pi + 
pi)j1 and ESR frequency uj. We take T = 70 mK, Ap/e = 
6 pV, B z = 0.5 T, g = 2, T x = 1 /is, T 2 = 100 ns, 71 = 
5xl0 6 s"\ and 72 = 571, i.e., A z = lOfcT and A/i = kT. The 
width of the sequential tunneling peaks in I(p) is determined 
by the temperature, see Eq. (|3lh . (a) The current J(/i, uj) 
shows a spin satellite peak near /1 = Es~ A z (for E-\ =0) due 
to the ESR field. Note that the spin satellite peak is slightly 
shifted from this position [see Eq. (p5[)], which is indicated 
by the line at Es — A z (light gray line) in (a). Here, B° = 
1.4 G, i.e., W™ ax = 71 at resonance and p = A Si . (b) The 
current for W u = (dotted), 71 /5 (solid), 71 (dashed), 
971 (dash-dotted). The position of the spin satellite peak as 
function of Wui is shown as black dots and the connecting 
solid line. 



can be measured. Such a measurement could provide a 
useful technique to study g factor modulated materials, 
where the g factor can be controlled by shifting the equi- 
librium position of the electronsJn the dot from one layer 
to another by electrical gatingO Note that measurement 
of the peak position would also allow to access the Stark 
shifts [Eq. (§)]. 

We consider now the analytic expression for the current 
J, as given in Eq. (Al), for the regime of the spin satellite 
peak. In this regime, Agf — p\ = Asi + A z — pi — 
Ap > A z - Ap M A z > kT, and thus //(A ST ) = 0, 
W s * — 0, and Wl s = 7j . For simplicity, we consider 
7; = 7/ = 7j here (cf. Sec. |Vl] for pumping due to 7/ ^ 
7/ ) . The expression for the stationary current [Eq. ([Al])] 
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considerably simplifies to 

I(cj, (i) =2e(W iT + W w )7i7 2 [/i(A si )-/ 2 (A S ;)] 

x{ (2 7 - W n - W u ) [71/1 (A Si ) + 72 / 2 (A Si )] 

+4 7 (W T1 + W iT +2W W )} . (29) 

For a plot of / vs u> anclu and some explanations of its 



characteristics, see Fig. 



B. Spin Decoherence Time T2 

Around the spin satellite peak, it is possible to mea- 
sure W u via the current and thereby access the spin de- 
coherence time of the spin-i state on the dot. For this, 
we identify a regime where the Rabi spin-flips on the 
dot become the bottleneck for electron transport through 
the quantum dot such that the current becomes pro- 
portional to the spin- flip rate W u - For kT < Afi and 
W™ &x < max{Wfi, 71} we obtain for the stationary cur- 
rent [Eq. ph], 



JH = 



2e7i7 2 (V^ iT + Wj) 
71(71+72) + (7i + 2 72) ' 



(30) 



see Fig. [| We have used W^ < W^i here. In the 
linear response regime, kT > A/i, and for W™ ax < 
max{W-fj,, 7/1 (Asi + A/x/2)}, the current is 

= 'WiWu+WjAy cQsh _ 2 ( A sl fi 



2{ 11 + l2 )kTh{T) 



2kT 



(31) 

The current I(fi) shows the standard sequential tunneling 
peak shape, determined by the usual cosh dependence on 
temperature, which is slightly modified by 

h(T) = 2W n + (2 7 -W n ) A(A Si + A^/2). (32) 

Most importantly, the current I(u>) of the satellite peak 
[Eqs. ( J30f ) and (|3l])] is proportional to the spin flip rate 
Wu. Thus, or equivalently I(B Z ), have a Lorentzian 
shape with resonance peak at lu = A z of width 2 Vj_-|- . 
Since Vjj > 1/T 2 , this width provides a lower bound on 
the intrinsic spin decoherence time T 2 of a single dot- 
spin. For weak tunneling, 71 < 2/T 2 , this bound satu- 
rates, i.e., the width 2V^ becomes 2/T 2 . Note that also 
Eq. ( |44| ) (see below) shows resonant behavior, i.e., a lower 
bound for T 2 can also be measured via a current due to 
pumping. 

We point out the similarity of our proposal to ESR 
spectroscopy]!!] where absorption or emission linewidths 
of the ESR field provide information on decoherence. 
In contrast to these techniques, we are considering here 
linewidths in resonances of the current, which allows us 
to access even single spins, since very low currents can 
be measured accurately. 

For Eqs. (^0|) and (|3l]) we have assumed that W u is 
small compared to the tunneling or the spin relaxation 



0.1 



< 



2/T 2 



14.75 A z /2n 14.76 



u/2tt 
[GHz] 



FIG. 3: The stationary current 7(w) [Eq. (gC()] for kT < 
Au, B z = 0.5 T, B° = 0.45 G, Ti = 1 us, T 2 = 100 ns, 
7i = 5 x 10 6 s" 1 , and 72 = 571, i.e., satisfying W" lax < 
71 < I/T2. Here, the linewidth gives a lower bound for the 
intrinsic spin decoherence time T2 (shown schematically by 
the arrow), while it becomes equal to 2/T2 for B® = 0.08G and 
WC ax < 7i = 5 x lO'^" 1 < 2/Ta, where I(lo = A*) w 1.5fA. 



rates. Therefore, we have neglected the contributions of 
W u in the denominator of these expressions. To take 
these contributions into account, we note that W u /{a + 
W u ) as a function of to is still a Lorentzian, but with an 



increased width w = 2V^-|-yl + W™ '- 
current I(u>) has the linewidth 



/a. Therefore, the 



w = 2V in l 



W™ ax (37i+4 72 ) 



71(71+72) + Wfx(7i + 2 72) 
for kT < Afi [Eq. ©], and 



(33) 



w = 2V LU /l + W l 



[A-h{A Si +A^2)]/h{T), (34) 



for kT > Afi [Eq. @]. Since the linewidth is increased 
by this correction, the inverse linewidth is still a lower 
bound for T 2 . 



C. Universal Conductance Ratio 

For increasing W^, the satellite peak in the current 
J(/x) increases while the main peak decreases, as shown 
in Fig. §(b). Further, as function of kT, the peak is 
slightly shifted. Explicitly, for jj = 7, , and A z > Afi, 
kT, we find from Eq. ( p9| ) the position of the satellite 
peak 

kT f iy n /2 + ^ iT +3^/2 + 7 

flESR = A Si - —\0g< — — — — 

2 I W n + W ir + 2W W 



The position of the main peak is 



kT f W n + 2W U + 3W U + 2 7 
Ho = A 5T + — log I w w 2Wu+2 ^ 



(35) 



(36) 



An experimentally accessible quantity is the ratio of the 
two current peaks or, equivalently (for linear response 
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0.71 



r(wg 



WJ 7 

FIG. 4: The current ratio r of the main and the satellite peak 
as a function of the effective spin- flip rate W u [Eq. @]. The 
dashed line shows the saturation of r for Win > 7 at the 
universal conductance ratio ro » 0.71 [Eq. (pq)]. 



Ap < fcT), the ratio of the conductances r(W 7 ^) = 
I(hesr)/I(po) = G(hesb)/G(p ). For this, we evalu- 
ate the stationary current at the gate voltages defined 
by Eqs. @ and @), and find, for A/j < kT and 



and p T = p| = 77 Wout, where 77 = l/(Wi n + 2W out ) is a 
normalization factor, ?7esr at the satellite peak and 770 at 
the main peak. The rates Wi n ( ou t) include all processes 
of electrons tunneling into (out of) the dot. Note that at 
the satellite peak, fj, = /jesr, a spin- up electron tunneling 
from lead 1 is the only process where an electron tunnels 
onto the dot, i.e., W- m (fiESR) = 7i, whereas at the main 
peak, fx = fio, the only tunnel process out of the dot 
is an electron with spin down into the right lead, i.e., 
W ou t(no) = 72- At the satellite peak, both spin up and 
down electrons can tunnel from the dot to lead 2, thus 
the current is given by /(/Zesr) = 2j2Ps = 27i7 2 7?esr, 
with 77ESR = l/(37i + 472). At the main peak, electrons 
can tunnel from lead 1 onto the dot, and the current is 
^(Mo) = 7i (PT + Pi) = 2 7i72%, with 7/0 = l/(27i + 372). 
Thus, the conductance ratio is given as r$ = t/esr/^Oj 
and we immediately obtain ro(oo) in accordance with 
Eq. (|3^). Therefore, the reason for ro < 1 is that at 
the satellite peak three out of four tunnel processes con- 
tribute to Wout, and thus t/esr < 770, while only one 
contributes at the main peak. 



2W U 



M' , 



21V„+27 



(37) 



IV. EVEN-TO-ODD SEQUENTIAL TUNNELING 



see Fig. [|. On the one hand, for small spin-flip rates, 
W u < 7, the ratio r is 4W u /7, i.e., at ESR resonance 
r(B°) = (gp B Bl) 2 /{2V ir y). If the tunneling rates and 
field strengths are known, this provides a further method 
for measuring a lower bound of the single spin decoher- 
ence time. On the other hand, this peak ratio [Eq. (37)] 
can be used to measure the ratio W u /j, useful for esti- 
mating the additional peak broadeni ng due to other lim- 
iting processes, as discussed in Sec. [II B , cf. Eqs. ( |3^ ) 
and @. 

It is noteworthy that this ratio saturates for W u ^> 7 
at the universal conductance ratio 



ro = 



2y/6 



7 + 4V3 



0.71 



(38) 



For a larger bias, but still Ap < A z , and for W u ^S> 7, 
the ratio becomes 



r ° (It) = 



V3 + V2e^ 



71 



U/2 + VSe^r 



72 



^2 + \/3e2fcT^ 71 + ^V3 + 2 e2ir^ 72 

(39) 

For 71 =72, the numerical value of ro remains 0.71 for 
all values A//. Generally, ro is between 2/3 (for 71 3> 72) 
and 3/4 (for 71 <C 72), where ro takes these extremal 
values for A/i > kT. 

Note that the current at the satellite peak is never 
larger than at the main peak. This asymmetry is best 
explained in the limit A/i > kT, when the ratio becomes 
r (oo) = (27 1 + 37 2 )/(37i+47 2 ). Since W w > 7, the Rabi 
spin flips equilibrate the populations p\ and p[. Thus, 
the stationary populations of the states are ps = 77 Wi n , 



Up to now we have considered sequential tunneling cur- 
rents with odd-to-even transitions of the number of elec- 
trons on the dot. Now we consider a different filling on 
the dot, with even-to-odd transitions. The state with N 
even is \S) (involving different orbital states as for \S)), 
and the states with N+l are | f ) and | J,). This system can 
be described with the same formalism as before, but with 
the tunneling rates — J2i W§ , , W^§ — J2i W\g ' 



w sl = i}[i-MA lS )] 



'I "si 
w ls 



7; /z(A 



is) > 



(40) 



and with W§^, W^§, Wg^, and W^ s defined analogously. 
The master equation of this system is given by Eqs. (|l6|)- 
( f2l| ) upon replacing the subscripts S by S. Since W\ § de- 
scribes an electron tunneling onto the dot, whereas Wis 
describes an electron tunneling out of the dot, the sta- 
tionary current through the dot is given by Eq. (|2^) after 
changing its sign and replacing the subscripts, resulting 



h = -e (W* s + W 2 lS ) Ps + e W| T p T + e W| ; p 4 . (41) 



By comparing Eqs. @, © with @, and Eqs. (|22 
with Q, we find that the formulas for the current are 
modified by the replacements fi(Asi) — ► [1 — fi(Aig)], 

7/ — > 7j , ij —> — Jj , and analogously for opposite spins. 
For completeness, we give in Appendix the formula for 
the stationary current 1% [Eq. (A3)], which is obtained 



by applyi ng the above replacements to Eq. (Al). 

In Sec. Ill B] we have identified the regime of the spin 
satellite peak, which can be used to measure the deco- 
herence time T2. For the setup considered here, an anal- 
ogous regime is p\ > Ajj > P2 > Ajg, see Fig. ||(a). 
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FIG. 5: (a) Setup for measuring T2, with > E± > [12 (for 
Eg =0). A lower lying state occupied by a singlet (corre- 
sponding to state \S)), illustrates the antiferromagnetic fill- 
ing of the dot. (b) Dot which should act as spin filter, al- 
lowing only spin f to pass. However, in the setup (b), the 
singlet-triplet spacing Ef — Eg is too small compared to 
Afi = fix — /J.2- Here, if the initial dot state is |f) (shown 
in gray), an electron with spin J. from a lower- lying state can 
tunnel onto the right dot, leaving a triplet on the dot (black), 
thus the spin filter does not operate properly. This problem 
disappears if the number of electrons on the dot can be re- 
duced down to zero. 




FIG. 6: Spin inverter setup, where the ESR field generates 
spin flips on dot 1, and the (additional) dot 2 acts as a spin 
filter, allowing only spin J, electrons to tunnel into lead 2. We 
consider the regime |£dd| < |toL 2 |, E\ ~ E%, A\ 96 A\, and 
Eg > m > E' s - A\, for % = 1, 2. The allowed transition 

sequence is schematically given by t©© 2 ^ t©© 2 

-^@© 2 ^©!@2 ^©©2 1 (see text), where "W 
means a coherent tunneling process. 



The current at the spin satellite peak is then given by 
Eqs. © and @ in the corresponding regimes, after 
interchanging 71 with 72, replacing fx — ► (1 — and 

For antiferromagnetic filling of the dot, one can use 
particle-hole symmetry to show that the two cases, odd- 
to-even and even-to-odd transitions, are equivalent. In- 
deed, the tunneling from, say, a spin ] electron from the 
dot into the lead, ||) — * \S), can be regarded as a spin 
f hole which tunnels from the lead onto the dot, which 
was initially occupied by a spin J, hole and now forms a 
hole singlet, i.e., | ih) — > \Sh). With this picture in mind, 
above modifications become obvious. 



V. SPIN INVERTER 



In this section we describe a setup with which spin- 
dependent tunneling, 7^ 7^ jj, can be achieved. Alter- 



natively, spin-polarized leads (see Sec. IX for details) or 
spin-dependent tunneling barriers could be used. This 
setup, shown in Fig. ||, consist of two dots, "dot 1" and 
"dot 2" , which are coupled in series with inter-jdpt tun- 
neling amplitude £dd ■ Dot 2 acts as a spin filtero and is 
coupled to the lead 2 with tunneling amplitude £dl 2 ■ We 
write the Zeeman splitting A d z , the energy E% of state 
\n), and the chemical potential Ag a with an index for 
dot d — 1, 2. We assume that dot 2 remains unaffected 
by the ESR field, which can be achieved e.g., by applying 
B x and/or B z locally or with different g factors for dot 
1 and dot 2. This assumption is taken into account by 
choosing A J 56 A 2 Z . 



A. Spin Filter 

We briefly jceview the concept of using a quantum dot 
as spin filter,E-3 as it is important for the description of 
the spin inverter. If the dot is initially in state 1 1), only 
a spin I electron can tunnel onto the dot, forming a sin- 
glet. Most importantly, the Zeeman splitting in the dot 
should be such that A z > A.$i — (J.2- This ensures proper 
operation of the spin filter: because of energy conserva- 
tion only the electron with spin J, can tunnel from the 
dot to the lead, leaving the dot always in state | |) after 
an electron has passed. Therefore, the sequential tunnel- 
ing current is spin J, polarized. There is a small spin-f 
cotunpcling current, however, which is suppressed by a 
factoiEil 7max{fcT, Afi}/(E T+ -E s ) 2 . Note that for effi- 
cient spin filtering, it is favorable to have the singlet state 
|5) as ground state with an even number of electrons on 
the dot, since the denominator of the suppression factor 
can become large, i.e., Et + — Es > A z . Otherwise, if 
the triplet state \T + ) — \ ft) is the ground state, only 
spin-| sequential tunneling current can flow through the 
dot. However, the spin- J, cotunncling current involves the 
triplet state |T ) = (|U) + I4T>)/V2, and the suppres- 
sion factor is given by 7max{fcT, A^i}/(A^) 2 , j-e., the 
cotunneling current is not suppressed efficiently.^ 



B. Implementation of Spin Inverter 

For implementations of the spin inverter, the Zeeman 
splitting in dot 2 should be such that A\ > A^ — (12, 
ensuring that dot 2 acts as a spin filter. The coupling 
of dot 2 to the lead shall be strong such that electrons 
escape rapidly from dot 2 into lead 2. This leads to res- 
onant tunneling with resonance width T2 = 27WjJ£dl 2 | 2 - 
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We require T 2 < A 2 



si 



- /i 2 , i.e., that the broadened level 
This excludes contributions from 



of dot 2 is above /z 2 
electrons tunneling from lead 2 onto dot 2, as shown in 
Ref. ||. 

We calculate the rates 7^ and 7^ for tunneling top, 
dot 1 via dot 2 into lead 2 in a T-matrix approach. L3l£| 
We use the tunnel Hamiltonian Ht — i?DD + Hdl 2 > 
where -Hdd describes tunneling from dot 1 to dot 2 and 
i?DL 2 from dot 2 to lead 2. The transition rates are 

Wfi — 27r|(/|T(£j)|i)| 5(ef—Si), where lead 2 is initially 
at equilibrium and with the T matrix 



T( £i )= Urn H T J2 

n=Q 



IT] 



H, 



dot 



Hi 



~ I J J' 



lead 



(42) 

We take the leading order in i?DD and sum up the contri- 
butions from all orders in i?DL 2 ■ We then integrate over 
the final states in lead 2 and obtain the Breit-Wigner 
transition rate of an electron with spin | to tunnel from 
dot 1 to lead 2 via the resonant level Eg of dot 2, 



|i DD | 2 r 2 



si 



A| T ) 



(r 2 /2V 



(43) 



In the spin filter regime considered here, dot 2 is always 
in state |t)- Thus, tunneling of an electron with spin 
t would involve the triplet level Et + on dot 2, which is 
out of resonance, and thus 7^ is suppressed to zero (up to 
cotunneling contributions, see Sec. |5 F| ). The state of dot 
1 and the current through the setup is again described by 
the master equation [Eqs. (^6|)-((2J)] with the tunneling 
rates W§ 1 = Wh = = and Wh = 1 L - Thus, 



we can use all previous results for one dot in Sec. HI A, 
but with 7 2 4 -> 7- 1 , 7 2 T -> 0, and / 2 (A ST ) = 0. Note that 
even for zero bias Ap = 0, a pumping current flows from 
lead 1 via the dots 1 and 2 to lead 2, see Eq. (04) and 
Sec. jv|. We point out that this setup, see Fig. £3, acts 
as a spin inverter, i.e., only spin f electrons are taken 
as input (lead 1), while the output (lead 2) consists of 
spin I electrons. In particular, the spin inverter does not 
require a change in the direction of the external magnetic 
field. 



to obtain a directed current, the spin symmetry between 
lead 1 and 2 must be broken. This can be achieved by 
spin-dependent tunneling, 7/ ^ 7/, e.g., produced with a 
double-dot, see Sec. [v[ At the spin satellite peak and for 



zero bias, i.e., /1 =f-/ 2 , there is a finite current [Eq. ([A 2])] 
due to "pumping"c9 by the ESR source, 

I(w) = eiW^+W^M ~ 7i7 2 T )/i(A5i)[2 7 Vi(Asi) 



x (2 7 x - W n - WJ) + 4 7 (V^ Ti +W U +2WJ 



(44) 



Here, sgn(7|72 — 7i7 2 ) determines the direction of the 
current. Note that for spin-independent tunneling, 7, = 
7, , and the pumping current vanishes. 

VII. ROTATING ESR FIELDS 

It is interesting to study rotating magnetic fields in 
addition to linearly oscillating fields as studied above. 
With rotating fields, it is possible to calculate the time 
evolution of the density matrix of the dot exactly. In par- 
ticular, the stationary solution of the master equation is 
obtained in a controlled approach and no rotating wave 
approximation is necessary. However, rotating fields are 
experimentally more difficult to produce than linearly os- 
cillating fields. 

We consider a clockwise rotating field, described by 

#esr = -j A_l [<j x cos(wt) - o~y sin(urf)], (45) 

where A^ = 2g/iB-B^. Thus, for A^, = A^ we have 
chosen the amplitude of the rotating field to be only 
half the amplitude of the linearly oscillating field, since 
both lead to the same effective spin flip rate W u . Us- 
ing Eq. (|5|) we immediately obtain the master equation, 
which is given by Eqs. (|l6|)-([2l|) after the following re- 
placements. The last terms in Eqs. (n6[) and dl7|) become 



T-(A ± /2)Im[e 1 ' 
placed by 



"Pit 



respectively. Equation (|19|) is rc- 



VI. PUMPING 

The ESR field provides energy to the system by excit- 
ing the spin state on the dot. When the dot is initially 
in the excited state ||), a spin up electron can tunnel 
onto the dot, followed by the spin down tunneling out 
of the dot. In total, the Zeeman energy A z is gained. 
This energy input can be exploited to induce a current 
through the dot, even at zero bias A/j = 0. However, 



Pu = ~iA zPi] + i_ e -^ t ( PT - P[ ) - V u m . (46) 

We transform to the rotating frame, |f) = e 4 "*/ 2 |t), 
and ||) r = e-^*/ 2 |4), such that Pl1 = e'^pW. This 
transformation removes the time-dependence of the coef- 
ficients in the master equation, which we shall now write 
as = A4pv- The equations for psi and psi decouple 
and we write the remaining part of the superoperator M. 
as matrix in the basis {pp p^, p r s , Re[/5^], Im[p^] 
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M 





w n 







Wn 


-(W u +Wsi) 


W iS 







Wsi 


-(W T s + W;s) 














-v L1 


Ax/4 


-Ax/4 





-(A,- 



-Ax/2 \ 
Ax/2 1 




(47) 



The master equation can now be solved exactly by cal- 
culating the eigenvalues Xi of KA. Since the total prob- 
ability is conserved, J2 n Pn = = E™^™?™. where 
n is summed over the diagonal elements, and m over di- 
agonal and off-diagonal elements of po- By considering 
linearly independent initial conditions for po, we see that 
Y^ n -Mnm = 0, for every m. Thus, adding up the rows 
in Ad for the diagonal elements of pd gives zero, which 
is satisfied explicitly by adding the first three rows in 
Eq. (f47|). Therefore, M does not have full row rank and 
there is an eigenvalue Ao = with eigenspace describing 
the stationary solution. The eigenvalues of M. are 



o,-vit,-3w; 



^(E vv + nT±\/(^-V|T) 2 -Ai , ( , 



with Sh/ — W + Wfi + Wif , and where we have consid- 
ered W = Wsi = Ws± = Wis = Wis, and resonance 
A z = lu for simplicity. If all Xi are different, the tirpft 
evolution of the density matrix is pr>(t) = ^ Ci£ Xit Pi*^ 
The decay of the contribution of the eigenvectors pi is 
exponential and generally all decay rates A^ are involved. 
Further, we see from the last two eigenvalues in Eq. 
that the decay rates of pd may be a nontrivial function 
of the rates involved in the master equation. This should 
be kept in mind when one uses time dependent ensemble 
properties, i.e., po(t), to measure intrinsic rates, e.g., T± 
and Ti. We point out that the presence of very small 
decay rates does not necessarily prevent a decay of the 
initial conditions. If, say, the tunneling rates are smaller 
than the spin relaxation rate, W <C W\\ % it would be 
interesting to study a density matrix which is described 
as a linear combination of the eigenvector with eigen- 
value — 3W [Eq. (|4S|)1 and the stationary solution po, 
i.e., Pd(£) = Po + ce~ 3Wt P3W, which is independent of 
W-fj_. However, such an initial condition always contains 
contributions from state \S) such that, in particular, it 
is not possible to construct an initial spin-^ state which 
would decay only with the slow rate 3W. 

The (exact) stationary solution of the master equation 
can be readily obtained from Eq. (|47|). By eliminating 
Pll from the coupled equations, we obtain the effective 
spin flip rate 



Ai 



IT 



8 (w — A z ) 2 + 



(49) 



which is equivalent to Eq. ([25|). Thus, all the results for 
the stationary currents from Sec. Ill apply and are exact 
for the case of rotating magnetic fields. 



VIII. COTUNNELING 



We now consider the cotunneling regime 



A 



ST i 

Ag i > Mi, fJ-2 3* Ei, Ei, where the number of electrons 
on the dot is odd, thus the state on the dot is described 
by |t) and \\). The leading order tunnel processes is 
now the tunneling of electrons from lead I onto the dot, 
forming a virtual state |n), followed by tunneling into 
lead V . The spin state of the dot changes a — > a' . This 
(48) process is called elastic cotunneling for a — a' and in- 
elastic cotunneling for a ^ a' . Note that in the absence 
of an ESR field, the dot relaxes into its spin ground state 
and no inelastic cotunneling processes, exciting the dot 
spin, occur for A/x < A z . However, if an ESR field is 
present, the dot-spin can be excited by spin flips. Then, 
inelastic cotunneling processes, which relax the dot-spin, 
can occur. These processes either contribute to transport 
or produce a particle-hole excitation in lead 1 or 2 [see 
Fig. @(b) and (c)]. 

These cotunneling rates are calculated in a "golden 
rule" approach ,Ej which is. known to be consistent with 
a microscopic derivationJU 



W: 



VI 



2ttp 2 defi(e) [1 - //(e- A ctV )] 



E 



t t* 



(50) 

where the possible spin-dependence of v has been ab- 
sorbed into t, and A a > a = E a i — E a is the change of Zee- 
man energy on the dot, and A na = E n — E a is the energy 
cost of the virtual intermediate state. Here, ti an are the 
tunneling amplitudes, where tiis — tl has already been 



introduced in Eq. (lla). The cotunneling current through 
the dot can be calculated by summing up the contribut- 
ing tunneling rates, as we have done for Eq. (E2), 



^ct — 



e£(Wa 



(51) 



We point out that by treating the cotunneling processes 
with golden rule rates, only classically allowed dot-states 
are considered. Thus, the number of charges on the dot 
is fixed and no charge can temporarily accumulate as for 
sequential tunneling. In particular, we have neglected 
quantum charge fluctuations on the dot. Therefore, 
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FIG. 7: Cotunneling processes involving \S) for A z > A/i. 

(a) Elastic cotunneling. The cotunneling sequence 1(T) — » 

(fi) — > (T)l , involving the virtual state IS 1 ) on the dot with 
virtual energy cost Ag-j — /ii. An equivalent process is possi- 
ble when the initial and final dot state is ||), however, with 
a virtual energy cost reduced by A z . These elastic cotunnel- 
ing processes contribute to transport and to spin decoherence, 
while they do not contribute to spin relaxation (i.e., Ti). (b) 
Inelastic cotunneling from lead 1 into lead 2 via the sequence 

T(D -> (Q) -> © I • Note that tunneling of an electron from 
lead 2 into lead 1 is also possible, since the energy gain A z 
from the dot relaxation is larger than the bias A/Lt. (c) Inelas- 
tic cotunneling, where only one lead is involved. The process 
shown here leads to a particle-hole excitation in lead 1. While 
it does not directly contribute to transport, it contributes to 
spin relaxation and spin decoherence of the dot. 



PIT = -i^zPi] +i^cos(W)0 T 



pi)-vR t piv (55) 



Note that away from the sequential tunneling regime, the 
master equation becomes much simpler while the formu- 
las for the rates are more involved. 

For the time-averaged current we evaluate the station- 
ary solution of the master equation in the rotating wave 
approximation (see Sec. Ill) for linearly or exactly (see 
Sec. VII) for circularly polarized ESR fields. This yields 

Eqs. (|5j) and (fi|), respec- 
We obtain 



an effective spin-flip rate W L _ 
tively] and eliminates Eq 



Pi = 



W w + W 



11 



2W UJ + W 



It 



w- 



n 



(56) 



TI 



and /Of = l—pi- We consider the case close to a sequential 
tunneling resonance (but still in the cotunneling regime) , 
As CT — pi < Et + — Es, such that the virtual energy cost 
of an intermediate triplet state is much higher than that 
for a singlet state. Since (Et + — E a — p)/[Es — E a — 
p) < 1, with fi = (/ii + /Lt2)/2, we have to consider only 
cotunneling processes involving state \S) in Eq. (|50|). For 
Ap, kT < As i —pi < Et + — Es, the relevant elastic rates 
are 



within our master equation approach for cotunneling, the 
charge currents in both leads are equal, I\(t) = hit)- 
This equality is valid for "coarse-grained" expectation 
values of the current (and other physical observables). 
In this approximation, one smoothens out the quantum 
fluctuations by averaging over the short-time behavior, 
i.e., one considers only the behavior on time scales larger 
than the lifetime l/(Ag| — p) of the virtual states on 
the dot. However, when the charge imbalance due to 
the virtual states is taken into account in a microscopic 
treatment, one can find pronounced peaks in the noise 
S(u>) for \lo\ corresponding to the virtual energy cost, as 
it was shown in Ref. ^5|. 

The inelastic cotunneling provides spin relaxation pro- 
cesses in addition to those contributing to T\, totaling 



^Ti + J2u> W-jv ^ or P rocesses w ith I' = I, 



particle-hole excitations are produced in lead I. We are 
interested in the regime A/i < A z , where (inelastic) co- 
tunneling does not excite the dot-spin, i.e., Wy? — W^. 
In analogy to Eq. ([l5]) , we take a phenomenological total 
spin decoherence rate 



t/CT 

TT 



1 



E 

11' aa' 



W: 



i'i 



(52) 



where all spin relaxation and tunneling processes are 
taken into account. The master equation for the dot in 
the cotunneling regime and in the presence of a linearly 
polarized ESR field becomes 



Pi = -wffn + wffpi 



P\ 



A x cos{ojt) Im[p; T ], (53) 
Wf? p T - Wff Pl + A x cos (tot) Im[ PiT ] , (54) 
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7i 72 



A/i 



2tt (Asa-p) 2 ' 
The inelastic rates are, for lead indices I, I' = 1, 2, 

A z + (l'-l)Ap 



7i 72 

2tt ~(A sl -p)(A sl + A z 
A, + (/'-/) A// 21 
A~p Wl1 ' 



P) 



(57) 

(58) 
(59) 



where Eq. ( |59| ) is valid for A z < A$i — p. Note that 
for An < A z the inelastic rates can be much larger (by 
a factor of A z /Ap) than the elastic ones, while their 
contribution to the current, — W^f — 2Wfl, is of 
the same order as for the elastic rates. 

For VF™ ax , Wi;i < W+i, we obtain the cotunneling cur- 
rent from Eqs. <M) and 



e A/X7172 
2^(A s ^p) 2 

e A^7r7 2 
2tt (A st - pf 



A 



si - P 



A 2 



Asi 



- P 
A/i 
2A7 



Ast - p 



(60) 
(61) 



The first term in Eq. (|6C|) results from elastic cotunnel- 
ing with spin ground state | f ) on the dot. The second 
term represents the increased current if the spin is flipped 
into state | j) before cotunneling occurs, since then both 
elastic and inelastic cotunneling processes contribute to 
the current. The current icT is proportional to W u , up 
to a constant background, and thus shows, as a func- 
tion of to, a resonant peak at w = A z of width 2Vj_-j- . 
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Thus, the intrinsic spin decoherence time Ti is accessi- 
ble in the cotunncling current as well as in the sequen- 
tial tunneling (see Sec. Ill B ) . Generally, the cotunneling 
current is much smaller than the sequential tunneling 
current, and thus it might seem more difficult to de- 
tect T2 in the cotunncling regime. However, since cur- 
rent and decoherence rate due to tunneling are propor- 
tional to 7 2 , the small currents can be compensated by 
choosing more transparent tunnel barriers, i.e., larger 7. 
Then, current and decoherence rate in the cotunncling 
regime can become comp arabl e to the sequential tunnel- 
ing values given in Sec. IIIB. For illustration we give 
the following estimates. For B z = 1 T, 5° = 2 G, 
g = 2, 71 = 72 = 5 x 10 9 s _1 , Ti = 1 (is, T 2 = 100 ns, 
Asi — P — and A/i = A z /5, the cotunneling cur- 
rent as function of the ESR frequency u> is 0.17 pA 
away from resonance and exhibits a resonance peak of 
/max = o 31 pA; with na if_ wi dth = 3.41 x iq7 s -l. 



IX. SPIN READ-OUT WITH SPIN-POLARIZED 
LEADS 

An electron spin on a quantum dot can be used as a 
single spin memory (or as a quantum bit for quantum 
computationEil) , if the spin state of the quantum dot can 
be measured. It was shown that a quantum dot con- 
nected to fully spin-polarized leads, A^ oads > ep > A z , 
can be used for reading the_spin state of the quantum 
dot via the charge current E3 Such a situation can be 
realized with magnetiCpjsemiconductors (with effective g- 
factors exceeding 100)0 or in the quantum Hall regime 
where spin-polarized edge states are coupled to a quan- 
tum dot £3 If the spin polarization in both leads is j, no 
electron with spin J. can be provided or taken by the leads 
(since = 0), and the rates Ws\ and W^s vanish. Thus, 
if the dot is initially in state 1 1), no electron can tunnel 
onto the dot (the formation of the triplet is forbidden 
by energy conservation) and I = 0, up to negligible co- 
tunncling contributions. However, if the dot is in state 
1 1), a current can flow via the sequential tunneling tran- 
sitions TQ ->■ (Q) -> Qt . Therefore, the initial spin 
state of the quantum dot can be detected by measuring 
the current through the dot. Note that for this read-out 
scheme, it is not necessary to have A z > kT on the dot, 
the constraint of having spin-polarized leads is already 
sufficiently strong. 

In the stationary regime and for A 2 > kT, the current 
becomes blocked due to spin relaxation (Wjj,). However, 
this blocking can be removed by the ESR field producing 
spin flips on the dot (with rate W u ). For W u < Wn, 
this competition leads again to a stationary current with 
resonant structure, 



I(w) = e{Wu+W u ) 



7i72 



72 W n + (71 +72) Wi T ' 



(62) 



from which Vjj (and 1/T 2 ) can be measured. Note that 
the relaxation rate W-\[ is rather small, thus only small 



ESR fields can be used, which leads to small currents. 



A. Counting Statistics and Signal-to-Noise Ratio 

We analyze now the time-dynamics of the read-out of 
a dot-spin via spin-polarized currents. The goal is to 
obtain the full counting statistics and to characterize a 
measurement time t mcas for the spin read-out. While we 
have considered only averaged currents so far, we now 
need to keep track of the number of electrons q which 
have accumulated in lead 2 since t = O.cZI The time evo- 
lution of pp(g, t), now charge-dependent, is described 
by Eqs. (|Iq)-(pi|), but with replacements W 2 S ps(l) — * 



W? q ps{g-l) in Eq. @, and W? 



SiPl(l) " Si 1 
we consider the distribution function 



Wl Pl (q+1) 



is 

in Eq. (|18|). Next, 
Pi(q, t) — YlnPnili t) that Q charges have accumulated 
in lead 2 after time t when the dot was in state \i) at t = 0. 
For a meaningful measurement of the dot-spin, the spin 
flip times W^, W^ 1 , and 1/A X must be smaller than 
imeas an d are neglected. Eqs. (|l6|) (pT ) then decouple ex- 
cept Eqs. (0) and (H), which we solve for — 1, and 
for pi = 1 at t = 0. The general solution follows by linear 
combination. First, if the dot is initially in state ||), no 
charges tunnel through the dot, and thus P\(q, t) = 5 q Q. 
Second, for the initial state ||), we consider kT < A/i 
and equal rates Wg, = W? s — W . We relabel the den- 
sity matrix p| (<?) —> Pm=2q, and p s (q) -> p m =2 q +i 1 and 
Eqs. ( |Tt| ) and (|l^) become 



pm = W{p m -i - p m ), 



(63) 



with solution p m (t) = (Wt) m e wt /ml (Poissonian dis- 
tribution). We obtain the counting statistics 



n . , (Wt) 2 "e- Wt f Wt 

P\ (q, t) = - — 7—T-, — 1 h 



(64) 



Experimentally, Pi(q, t) can be determined by time series 
measure ment s or by using an array of independent dots 
(see Sec. XA). The-Lrpprse signal-to-noise ratio is defined 
as the Fano factor ,OC?J which we calculate as 



Flit) 



(Wt) 2 ) 

(<?(*)} 



3-2e- 2Wt {4Wt+l)-e- 



4(2Wt-l + e- 



2Wt\ 



with Fi decreasing monotonically from ^(0) 



AWt 

(65) 
1 to 

Fi(t —> 00) = 5. Note that for dot-spin ||\_only weak 
cotunneling occurs with Fano factor F| = 1.E3 

If we are interested in the current and noise for l ong 
times t > W~ l , we can follow the steps used in Ref. |5(J. 
We decouple the differential equations with respect to 
q by taking the inverse Fourier transform, pT>[k) — 
J2 q e~ lkq pY)(q)- Note that, for k = 0, we recover the 
density matrix p^ = po(k = 0), where the accumulated 
charge is not taken into account. The probability Pi (q, t) 
is approximated by a Gaussian wave packet in g-space 
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with group velocity I/e 
width y/2F{I/e% and 



Wis), 



and 



F = 



(w sl f + (wf s f 



si 



(66) 



is the Fano factor .E3 However, within this approximation, 
valid for Wt > 1, we cannot access the short time behav- 
ior where only a few electrons have tunneled through the 
dot, which is of importance for the read-out process con- 
sidered here. 



B. Measurement Time 

Using the counting statistics, we can now quantify the 
measurement efficiency. If, after time t meas , some charges 
q > have tunneled through the dot, the initial state of 
the dot was | J.) with probability 1 (assuming that single 
charges can be detected via an SETej). However, if no 
charges were detected (q = 0), the initial state of the spin 
memory was | j) with probability 



1-P X (0, t) = 1 



-w{ s t 



W 



is* 



Wf s 



(67) 



which reduces to 1 — e~ wt (l + Wt), for equal rates. Thus, 
roughly speaking, we find that t me&s ^ 2VF -1 , as ex- 
pected, while the Fano factor is 0.5 < Fj, < 0.72. If, more 
generally, the threshold for detection is at to charges, 
to > 1, Eq. @ is replaced by 1 - Y^=o Pl(Q, t). _ 

We insert now realistic numbers to obtain an estimate 
of the fastest possible measurement time which can be 
achieved with this set-up. For a fast spin read-out, the 
tunneling rates and the current trough the dot should be 
large, limited by the fact that the conductance of the dot 
should not exceed the single-channel conductance e 2 /h. 
In the linear response regime and for a small bias A/i/e, 
the current is I = e^Ap/8kT < (Ap/e) x (e 2 /h), for 
jI =72- Thus, the tunneling rates are limited by 7^ < 
8kT/h = 1.76 x 10 11 (T/K) s _1 . For W = 7 T = 1.25 x 
10 10 s _1 (corresponding to kT < Ap, and a current / = 
1 nA), and m = 1, the spin state can be determined with 
more than 95% probability for a measurement time of 
^mcas = 400 ps, rand with more than 99.99% probability 
for i mcas = 1 nsO 



X. RABI OSCILLATIONS OF A SINGLE SPIN 
IN THE TIME DOMAIN 

A. Observing Rabi Oscillations via Current 

The ESR field generates coherent Rabi oscillations of 
the dot spin, leading to oscillations in pr>(t). Since the 
time-dependent currents I{t) in the leads are given by 
the populations p n (t) [Eq. (f22[)], current measurements 
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FIG. 8: Rabi oscillations of the electron spin on the dot in 
the time domain. We consider the regime at the spin satellite 
peak, As T > /ii > A Si > /12, (see Fig. |l]) and take Ti = 1 /is, 
T 2 = 300 ns, A x = 5W s .i (corresponding to B° = 10 G for 
g — 2), and = 1 at t = 0. During the time span shown here, 
less than 3 electrons have tunneled through the dot on aver- 
age. Here, the spin decoherence is dominated by the tunneling 
process, i.e., Wsi S> I/T2. (a) Spin-polarized leads with the 
only non-vanishing tunnel rates Wsi = Wis — 4 x 10 7 s _1 . 
The Rabi oscillations show up in (dotted), pi (dashed) 
and ps (full line), which is directly visible in the current, since 
ll(t) oc pi and/jW K Ps, for kT < A/i. In the inset, we show 
the case of large tunneling, Wsi ~ W| g = 1 9 s~ 1 S> A x . Asa 
consequence of the Zeno effect (see Sec. |X C| ) , the Rabi oscilla- 
tions are suppressed. Further, pi and ps are indistinguishable 
since | J,) and |S) equilibrate rapidly due to the increased tun- 
neling, (b) The time-dependent currents in unpolarized leads, 
h(t) = e r yi(pi — ps) and h{t) = 2e r y2ps, for kT < A/i, and 
jj = 7 '=4x 10 7 s" 1 , for 1 = 1,2. 



give access to these Rabi oscillations. First, we consider 
a dot coupled to unpolarized leads in the re gime of the 
spin satellite peak (see Fig. |l| and Sec. [II A ). For kT < 

A/i, the current in lead 2 is l2{t) = + I2) Ps(£h 
i.e., ps is directly accessible via measurement of 7 2 (i)-E3 
Further, for 7J = 7^, the current in lead 1 is h(t) — 
e li(Pl ~ Ps)i which gives access to Pi(t), if the ratio 
71/72 is known. We calculate the oscillations of Ii^it) 
explicitly by numerical integration of the master equation 
[Eqs. ©-©], see Fig. |(b). 

The measurement of p^> can be refined by using the 
spin read-out setup with spin-polarized leads (Sec. |lX|). 



For kT < A/i, the current is h(t) = l[ (t) = ejiqi(t) in 

lead 1, and I 2 (t) = J.f(i) = ej 1 2 p s (t) in lead 2.E3 Thus, 
the time-dependence of pi and ps (and also of = 1 — 
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Pl — ps) can be directly measured via the currents Ix,2, 
see Fig. |(a). 

Note that the electrons which tunnel ont o th e dot de- 
cohere the spin state of the dot (see Sec. |IIE| ). Thus, 
to observe Rabi oscillations in h,2(t) experimentally, the 
Rabi frequency A x must be larger than the coupling to 
the leads Ws\., otherwise the strong decoherence (equiv- 
alent to a continuous measureme nt) su ppresses the Rabi 
oscillations (Zeno effect, see Sec. X C ). Then, however, 
only very few electrons tunnel per Rabi oscillation period 
through the dot. To overcome the limitations of such a 
weak current signal and to obtain I\ t 2(t) experimentally, 
an ensemble average is required. 

There are two possibilities to obtain averages, namely 
using many dots or performing a time series measure- 
ment. First, many independent dots can be measured 
simultaneously by arranging the dots in parallel to in- 
crease the total current. For example, an array (ensem- 
ble) of dots and leads could be produced with standard 
techniques for defining nanostructures, or self-assembled 
or chemically synthesized dots could be placed within an 
insulating barrier between two electrodes. Second, time 
series measurement over a single dot can be performed. 
For this, the procedure of preparing the dot to the de- 
sired initial state, applying an ESR field and measuring 
the current has to be repeated many times (see Sec. [X A 
for counting statistics of the read-out process). Then, 
assuming ergodicity, the current average of all these in- 
dividual measurements corresponds to the ensemble av- 
eraged value. 



B. Decoherence in the Time Domain 

In Fig. ||, we plot the numerical solution of Eqs. (Jlq)- 
(pl|), showing the coherent oscillations of pd and for 
(a) spin-polarized and (b) unpolarized leads. The decay 
of these oscillations is dominated by the spin decoherence 
rate . Since this decay can be measured via the cur- 
rent, Pjf (and I/T2) can be accessed directly in the time 
domain (see also Sec. XI, Rcf. |53| and Fig. g). 



C. Zeno Effect 

When the rate for electrons tunneling onto the dot, 
Ws<t, is increased, the coherent oscillations of pj, pi be- 
come suppressed (see inset of Fig. ||(a)). This suppres- 
sion is caused by the increased spin decoherence rate 
[Eq. (|l5|)1 and can be interpreted as a continuous strong 
measurement of the dot-spin, performed by an increased 
number of charges tunneling onto the dot. This suppres^ 
sion of coherent oscillations is known as Zeno effect S3 
Since it is visible in pn, it can be observed via the cur- 
rents h,2(t)- 



XL PULSED ESR AND RABI OSCILLATIONS 

We now show that it is possible to observe the coher- 
ent Rabi oscillations of a single electron spin even with- 
out the requirement of measuring time-resolved currents. 
This can be achieved by applying ESR pulses of length t p 
and by measuring time-averaged currents (over arbitrar- 
ily long times). Then, the time-averaged current I(t p ) as 
function of t p gives access to the time evolution of the 
spin-Sjtate of the dot, for both, polarized and unpolarized 
leads. E3 In particular, since arbitrarily long times, and 
thus a large number of electrons, can be used to mea- 
sure /, the required experimental setups are significantly 
simpler compared to setups which aim at measuring time- 
dependent currents with high resolution. 

We assume a rectangular envelope for the ESR pulse 
with length t p and repetition time t r (thus t p < t r ). The 
time when no ESR field is present, t r — t p , should be 
long enough such that the dot can relax into its ground 
state 1 1), i.e., at the beginning of the next pulse we have 
Pi = 1. We calculate I{t v ) by numerical integration of 
the master equation [Eqs. (|K^-([l9|)] and by subsequently 
averaging the (time-dependent) current [Eq. (E2h] over 
the time interval [0, t r ]. The results are shown in 
Fig. |o|(b ) for u npolarized leads at the spin satellite peak 
(see Sec. |III A ), and in Fig. ||(c) for spii^polarized leads in 
the regime for spin-read out (see Sec. IX). In both cases, 
I(t p ) as function of pulse length t p shows the Rabi oscil- 
lations of the dot spin, i.e., the Rabi oscillations can be 
observed in the time domain even without time-resolved 
measurements. 

In addition to the exact numerical evaluation of the 
master equation (see Fig. |t]), we now give an approximate 
analytical expression for I(t p ). We first consider t he cas e 
of unpolarized leads at the spin satellite peak (Sec. Ill A ); 
for the case of spin-polarized leads see below. For this, 
we need to evaluate the time-average of Eq. (|2^) for kT < 
Ap, 



7 2 J — 



dt Ps (t) . 



(68) 



First, we consider times t with < t < t p , for which 
an ESR field is present, and pd oscillates with Rabi 
frequency A x (see Fig. ^|(a) for t < 200 ns). Quali- 
tatively speaking, when ps{t) is integrated in Eq. ( |6^ ) 
up to t p , the oscillating contribution averages nearly to 
zero, and we obtain a background contribution Iq ap- 
proximately proportional to + l^p/^r, i-e., linear 
in t p , in agreement with Fig. Mb). For experiments, 
this linearity of Jo provides a first check that t r is suf- 
ficiently long such that the dot has indeed relaxed into 
its ground state before the next pulse is applied. We 
also give an upper bound for Iq by using the inequality 
Ps < PT X = Wsi/{Wsi + Wsi + Wis)- This is seen 
as follows. For ps(t) > Pg 111 *, we would have ps(t) < 0, 
and thus ps{t') > p™ 3 *, for all < t' < t, which would 
be in contradiction to the initial condition ps(0) = 0, 
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FIG. 9: Single spin Rabi oscillations in the current I(t p ) 
generated by ESR pulses of length t p . Here, Ap > kT, Rabi 
frequency A x — 4 x 10 8 s" 1 (corresponding to g = 2 and 
B° = 20 G), 7i = 2 x 10 7 s"\ 72 = 571, 7\ = 1 /xs, and T 2 = 
150 ns. (a) Evolution of the density matrix for unpolarized 
leads where a pulse of length t p = 200 ns is switched on at t — 
0, obtained by numerical integration of the master equation 
[Eqs. (b) Time-averaged current I(t p ) (solid line) 

for unpolarized leads and a pulse repetition time t r — 500 ns. 
We also show the current where 71 and 72 are increased by a 
factor of 1.5 (dotted) and 2 (dash-dotted), (c) Time-averaged 
current I(t p ) (solid line) for spin-polarized leads, 7J = 2 x 
10 7 s" 1 , 7j = 57{, 7! 2 = 0. The pulse repetition time t r = 
10 ps is chosen larger than T%. Again, we show the current 
for tunneling rates 7J 2 increased by a factor of 1.5 (dotted) 
and 2 (dash-dotted). Note that in this figure t p < T\, i.e., 
for most electrons tunneling through the dot after the pulse 
is switched off, the linear background is negligibly small. 

hence indeed ps(t) < Ps ax - From Eq. (|6S|), we then ob- 
tain Jo < emin{7j, jI + "i\~\t p jt r . Note that for pulse 
lengths t p , over which the dot spin evolves coherently, 
ip7i ^ 1- Thus, by comparing the upper bound with 
Eq. ([39]) , we see that for 7] < 72 the background current 
Jo never becomes dominant. 

Second, we consider t p < t < t r , i.e., the ESR field 



is switched off, and the dot state relaxes into its ground 
state |t). Making the reasonable assumption that the 
tunnel processes dominate the spin relaxation, 7 > W-r;, 
we neglect W-\i here. We then calculate the contribution 
for t > t p to the integral in Eq. ( |68| ) analytically, and 
obtain (up to Jo) 

T i 

^l^rT^^W+^Wl a 1 -P^p)- (69) 

We now give a physical explanation for Eq. (^). We 
consider different tunneling events (after the pulse is 
switched off) and their contributions to the current, 
f* r dt ps{t). Since we assume that at t r the dot has 

relaxed into its ground state |t), and thus ps(t r ) = 
Pl(t r ) = 0, it is sufficient to consider only one pulse and 
to extend the upper integration limit to infinity. For 
the population p±(t p ) of state 1 J,), the only allowed tran- 
sition is \l) — > \S) (neglecting again the intrinsic spin 
relaxation rate Wfj,). Thus, eventually this population 
Pl will be transfered to ps and thus to the current. Note 
that sequences with \S) —> ||) contribute to the cur- 
rent at a later time again, since the only possible de- 
cay into the ground state ||) involves IS*). Therefore, 
concerning current contributions, we introduce the effec- 
tive population pj = P[ + ps, which is the probability 
that at some later time an electron can still tunnel from 
the dot to lead 2. This pi decays to state | ]) with the 
rate 7s = 7}" + 73, i.e., with the rate for the process 
\S) — ► |t)- In total, integrating over ps(t) for t > t p 
yieldsT dtp I (t p )e-~ 1st = [pi{t p )+ psM]/ls, and with 
Eq. ( pq ) we immediately recover Eq. (g9|), as expected. 

Next, we consider the case for spin-polarized leads. 
Here, no spin relaxation process due to tunneling occurs 
and the dot-spin can only relax via intrinsic spin flips, 
given by the rate W-fj, (corresponding to the relaxation 
time Ti; we neglect since -C Wn). Thus, we 
now consider the relaxation rate Wfj, instead of 7s. The 
relaxation occurs only from | J.) to 1 1), i.e., the roles of \S) 
and 1 1) are interchanged compared to the case for unpo- 
larized leads considered above. The above argument now 
applies analogously by considering the (spin-polarized) 
current in lead 1, jj(t) — p[{t). We obtain 

JT ^)«f w-ti-M^L (70) 

t r W^i 

with equality for t p <C T\. We point out that for 
7 1 ^> l/Ti, the decoherence of the dot-spin occurs much 
faster than its relaxation. Then, for pulse lengths t p , 
for which Rabi oscillations can be observed, are limited, 
^ Vj.f > 7i ^> Wfj.. In this case, the current con- 
tribution for t < t p can be neglected since they are sup- 
pressed by a factor of t p W^i -C 1 compared to the contri- 
bution for t > t p [Eq. @], see Fig. |(c). Note that for 
spin-polarized leads, the relaxation time, , is usually 
much longer than for unpolarized leads, Jg , thus the 
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required pulse repetition time t r > might become 

very long. However, if one chooses a pulse repetition time 
t r = e/7, for c > 1, and with the relevant relaxation rate 
7, the current is proportional to (l/t r ) J Q dte~ l1 = 1/c, 
i.e., independent of 7. Thus, roughly speaking, the slow 
relaxation rate in the case of spin-polarized leads has no 
influence on the attainable maximum current since the 
decay from ps and p^ is much slower and thus per pulse 
there are more electrons passing the dot. 

To conclude, we would like to emphasize again that the 
Rabi oscillations of the dot-spin can be observed directly 
in the time domain by using pulsed ESR and measuring 
time-averaged currents (see Fig. ^). Observing Rabi os- 
cillation s als o-allows to determine Ti in the time domain, 
see Sec. BcbI.LJ 



XII. STM TECHNIQUES AND ESR 

So far, we have considered a quantum dot coupled to 
leads. In this section, we would like to note that our 
description applies to more general structures showing 
Coulomb blockade-behavior, such as Au nanoparticlest3 
or Ceo molecules, e3 which has been observed with STM 
techniques. This justifies that instead of a quantum dot, 
we now consider a localized surface state or an atom, 
molecule, or nanoparticle adsorbed on a substrate. This 
particle can then be probed with the STM tip by measur- 
ing the tunnel current through the particle. The current 
arises from electrons tunneling from the STM tip onto 
the particle, and further tunneling, possibly through an 
insulating overlayer, into the bulk of the substrate. 

In standard STM theory, the tunnelins-irom the STM 
tip to the sample is treated pertubatively.t3 Evaluation of 
the golden rule matrix element, in the simplest model of a 
one-dimensional tunnel barrier, gives a tunneling ampli- 
tude, which is dominated by an exponential decay of the 
electronic wavefunction into the barrier, thus cx e~ Kd 
[cf. Eq. (fL3|)], with k = \/2m</>, tip-particle distance d, 
and barrier height <j> (roughly given by the work func- 
tion of the tip/sample). In particular, the perturbative 
description of STM is equivalent to our treatment of the 
tunneling Hamiltonian in first (sequential tunneling) or- 
der. Therefore, if the particle of interest shows Coulomb 
blockade behavior and has a spin-i ground state, the 
master equation [Eqs. (fifi|)-(|2l|)] applies. Thus, using an 
ESR field, coherent Rabi oscillations and the T2 time of 
the spin state of the particle can be accessed via the cur- 
rent. Further, if spin-polarized tips and/or substrates are 
available (spin-polarized STM) , such a particle can act as 
single spin memory with read-out via current. Note that 
the tunneling rates from the STM tip into the particle 
can be controlled by changing the distance d, thus the 
total decoherence Vj.f [Eq. (jig)], containing tunneling 
contributions, can be varied. This allows, e.g., to vary 
the current linewidth, 2Vj,j, (Sec. IIIB), and to suppress 



the Rabi spin flips for strong decoherence (Zeno effect, 
Sec. X C ) . One apparent restriction of atomic or molecu- 
lar systems is that it is difficult to apply a gate voltage to 
the particle, shifting its energy levels. However, the same 
effect can be achieved if the Fermi energies in the STM 
tip and the substrate can be shifted, such as by varying 
electron densities. 



XIII. DISCUSSIONS 

We have shown how single spin dynamics of quantum 
dots can be accessed by current measurements. We have 
derived and analyzed coupled master equations of a quan- 
tum dot, which is tunnel coupled to leads, in the presence 
of an ESR field. The current through the dot in the se- 
quential tunneling regime shows a new resonance peak 
(satellite peak) whose linewidth provides a lower bound 
on the single spin decoherence time T%. We have shown 
that also the cotunneling current has a resonant current 
contribution, giving access to Ti. The coherent Rabi os- 
cillations of the dot-spin can be observed by charge mea- 
surements, since they lead to oscillations in the time- 
dependent current and in the time-averaged current as 
function of ESR pulse length. We have shown how the 
ESR field can pump current through a dot at zero bias if 
spin dependent tunneling or a spin inverter is available. 
We have discussed the concept of measuring a single spin 
via charge in detail. We have identified the measurement 
time of the dot-spin via spin-polarized leads. Finally, 
we have noted that the concepts presented here are not 
only valid for quantum dots but also for "real" atoms or 
molecules if they are contacted with an STM tip. 
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APPENDIX A: STATIONARY CURRENT 

Here, we give the various formulas for the station- 
ary current through the dot in the sequential tunneling 
regime and in the presence of an ESR field. We have 
calculated the current by evaluating the stationary solu- 
tion of the master equation (Sec. Ill) and with Eq. (j2^). 
For odd-to-even sequential tunneling, the spin \ polar- 
ized current in lead 2 is 
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n = 



e 7 2 



(7 1 T E(- 1 )'^/KA5 i )/r(A ST )+E 
\ i,i' i 



2W U + W n + IF iT 



{(-l)' 7l T /KAs|) + 7 ; i [/ 2 (As i ) - /i(A ST )]} 



£ — 2 — {(-l)Si T //(A Si ) + 7 ; i [/2(A Si ) + /z(A 5T ) - 2/ 2 (A 5i )/ ; (A ST )]}J 
Y^lH {1 - [1 - //(Asi)] [1 - /i'(AsT)]} + £ + ^'-) K + 7f'[l - //(As.)]} 



(Al) 



The spi n { polarized current, I 2 , is obtained from 
Eq. (Al) by exchanging all T and J, in the numerator (the 
denominator remains unaffected by such an exchange). 

The currents in lead 1, are obtained from the for- 
t 1 

mulas for 1% by exchanging indices 1 and 2 and by a 
global change of sign. The charge current is I? 
and is equal in both leads, I = I\ = I2, due to charge 
conservation. For large Zeeman splitting A z > A/i, kT 
and aroun d the spin satellite peak, /ii > A$i > ^2 (see 
Sec. [II A ), we have fi(Asi) — 0, and the current is 



I = e (W u + W 1T ) [(7i T 7 2 T + 7i T 7 2 X ) h(^si) 
-(7i7 2 T +7i T 7 2 T )/2(A Si )" 



I, a^a' 

r 



i - w- 



U 



7 1 T /i(A Si )+7 2 T /2(Asi)' 



(A2) 
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+2 (W n + W u + 2W U ) ( 7 X + 7 T ) 



for which we have given special cases in Eqs. fl29j) , (|3C|), 
© and (@). 

For completeness, we also give the results for even-to- 
odd sequential tunneling, as discussed in Sec. [V. By 
applying the replacements given in Sec. IV. to Eq. (Al), 



we obtain the spin | polarized stationary current in lead 

2, 



In 



e7 i I .1 



(7i £(-l)SU(A iS )[l - /,,(A TS )] + E 2 ^ + ^ n + ^ iT {(-l)4/i(A;s) + 7, T [/a(A iS ) - /*(A T s)]} 
- E Wn ~ W " {(-lYrlMAs) - 7« T [/2(A iS ) + /;(A T s) - 2/ 2 (A ig )/ ; (A Tg )]}^ 
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plied to cancel the effect of dephasing. 

000 

At the end 

of the sequence, the dot spin (in the a z basis) can be mea- 
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sured by one of the techniques discussed in Sees 
and [xi|. Here, spin echo compensates for dephasing of an 
array of quantum dots in an inhomogeneous magnetic field 
B z , or in time series measurements (on a single dot), where 
the magnetic field B z fluctuates in time. More care must be 
taken for inhomogeneities in the g factor (Rabi frequency 
A x ), since then a pulse does not rotate the spin by an 
equal angle for every dot. To compensate this, one cpuld, 
e.g., use Carr-Purcell-Meiboom-Gill pulse sequences.u 
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